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1. INTRODUCTION 


A FEW years ago, Gell-Mann (1956) suggested a minimal principle for electro- 
magnetic interactions according to which the electromagnetic interactions 
between elementary particles can only be of the current-field type that occurs 
in quantum electrodynamics and cannot be purely of the magnetic moment 
type. In this paper we suggest two specific transformations, as additions 
to the now well-known mass reversal transformation, invariance under 
which will guarantee the minimal character of the electromagnetic inter- 
actions. Applications of these transformations to strong and weak interactions 
are considered. These considerations suggest a model of strong interactions 
and an isotopic spin scheme for leptons which are also discussed. 


2. THE Mass GAUGE TRANSFORMATION 


To understand the significance of this transformation, we notice that 
in quantum electrodynamics, the total mass of electrons minus the total mass 
of positrons is a constant of motion. This implies that the quantum electro- 
dynamical Lagrangian is invariant under the following gauge transformation 
of the electron field: 


te— e-tmey,, > he aap pee atte (1) 


where mg is the rest mass of the electron. If now we postulate that all the 
electromagnetic Lagrangians should be invariant under the transformations 


bh etMay,, ; ba > be" im, (2) 


where #, is the field and m, the rest mass of any fermion, we obtain the 
minimal principle of Gell-Mann apart from the case where a particle is 
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coupled to itself by a magnetic moment type of interaction. The conservation 
law implied here is that the total mass of fermions minus the total mass of 
anti-fermions should be a constant of motion in the electromagnetic inter- 
actions. In the Feynman picture, the principle may be stated as follows. 
As we follow a fermion line in a Feynman diagram there is nothing to prevent 
a coupling that will make a fermion become another fermion of different 
mass at a vertex, the only conservation law imposed at a vertex being on 
energy and momentum. However in quantum electrodynamics we notice 
that as we follow a fermion line, its mass is conserved. If we now insist that 
this be true in all electromagnetic interactions we obtain the result stated 
above. 


It is obvious that this invariance holds also for the strong pion-baryon 
interactions in the charge independence approximation where the mass differ- 
ences within a multiplet are neglected, except for the 2 — A — = coupling. 
In the global symmetry model [Gell-Mann (1957)] where the 2° — A® mass 
difference can also be neglected in the absence of K-meson interactions, this 
symmetry holds for all the 7-baryon couplings when the K-meson inter- 
actions are switched off. 


Turning now to the weak interactions, decays like p+-»e+ +7 +7 
show that the fermion mass conservation is violated in weak decays. To 
deduce the universal character of the weak interactions from the violation 
of the fermion mass conservation, we suggest the following argument. The 
a-baryon interactions, in contrast to the electromagnetic interactions, allow 
a change of the electric charge of the fermions. These couplings also exhibit 
a higher degree of symmetry in that the couplings are the same for the 
different charge states of the same field. Pursuing this analogy, we may 
expect the weak interactions to involve a universal Fermi coupling (as is 
experimentally true) since these interactions allow the fermions to change 
their mass states. In this picture, as far as the weak interactions are 
concerned, all the fermion fields are different modes of a single fundamental 
fermion field. 


As another application of the mass gauge transformation, we notice 
hat the K-meson interactions connect fermions of different masses. Hence 
we may similarly expect the K-mesons to be symmetrically coupled to all 
the baryons. This is the ‘cosmic symmetry’ model of strong interactions 
arrived at by Sakurai (1959) from a different point of view. Also, apart from 
the 2 — A—7a coupling, the 7-baryon couplings connect the same mass 
states of the fermions, neglecting of course the mass differences within a multi- 
plet. This suggests that perhaps the 2 — A — 7m coupling does not exist. 
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This has also been suggested by Feynman (1958). This will then ensure 
parity conservation for all the non-derivative pion interactions from CP 
invariance and charge independence in the conventional sense [Feinberg 
and Gursey (1959)]. The K-baryon interactions will also conserve parity 
because of their ‘ cosmic symmetry ’ [Sakurai (1959)}. 


3. THE Mass REVERSAL TRANSFORMATIONT 


Recently when studying the structure of Fermi interactions, Feynman 
and Gell-Mann (1958) showed that the Dirac equation can be written as a 
second order equation 


(ve sx, +) (35, — ™) ¥ =O. (3) 


If % is a solution of this new equation, %’ = y,% is also a solution. So is a 
linear combination of % and y,%. This is not true of the first order Dirac 
equation since y; does not commute with the Dirac Hamiltonian. ,% corres- 
ponds to the solution of the equation obtained by replacing m by —»m in the 
Dirac equation and this is equivalent to stating that % and y, are solutions 
of the Dirac equation of the second order. In considering the interactions 
of particles, Feynman and Gell-Mann suggested that these solutions of the 
second order equation should be treated as fundamental, of course imposing 
conditions depending upon the nature of the interactions. In the case of 
universal Fermi interactions, the linear combination of solutions used is 
4(1 + y;) # which are the eigen states of the chiral operator y+, | 
and Marshak (1958)]. 


Since % and y,% have opposite parities, we may insist that all the inter- 
action Lagrangians be invariant under the simultaneous reversal of the 
parities of all the fermion fields by the transformation 


pop = nysb (4) 


with |n|? = 1, since only the relative parities of fermions are physically sig- 
nificant [Wick, Wightman and Wigner (1950)]. This leaves the current-field 


type electromagnetic Lagrangian ieby WA, invariant provided that 
A,—A,’' = A,-. If we accept this transformation for A,, then the mag- 
netic moment type couplings are not. invariant under these transformations 
and are consequently forbidden. 





t This transformation has been discussed in various contents by Sakurai (1958), Okun (1958) 
and Salam (1959). We merely collect them here as a. single connected series of results. 
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-?>> . 
If we now require Yukawa couplings like igpy;r to be invariant under 


J 
this transformation, we find that ¢ must transform as 


oe > 


b> =—4. (5) 


With this transformation for $. it follows that there can be no derivative 
Yukawa couplings. Invariance under this transformation also forbids 
couplings between an odd number of spin zero bosons like A¢*. [see also 
Baym (1960)]. These results are true whether ¢ is scalar or pseudoscalar. 


4. THE tT, TRANSFORMATION 


The minimal principle for the electromagnetic interactions of baryons 
may be formulated by specifying the form of the corresponding interaction 
Lagrangians in the isotopic spin space. If py is the nucleon field, its electro- 


magnetic interaction can be written as ie/2 Yyyy(l +75)%yA,. Such an 
interaction is obviously not invariant under rotations about the first or second 
axis in the isotopic spin space and consequently it is not invariant under a 
reflection in the plane formed by the first and third axes of the isotopic spin 
space. Therefore if we define a parity for the fields for such a reflection, 
this parity is not conserved in electromagnetic interactions. This is ana- 
logous to what happens in weak interactions where space parity is not con- 
served. In the case of weak interactions, the universal V + A interaction 
can be arrived at by arguing that the interaction should be unchanged even 
if the fields that occur in it are taken over to a state of opposite parity indi- 
vidually, as parity is not conserved. [Sakurai (1958), Sudarshan and Mar- 
shak (1958)]. This is achieved by the transformation y—>’ = + y,y. 
In a similar manner, we can expect the electromagnetic interactions to remain 
invariant under a transformation which switches the parity of a field under 
the reflection in the isospin space referred to above to the opposite parity. 
It is important that this transformation can be applied separately to each 
of the fields in the interaction. Such a transformation for isospinor fields 
would be 


fb —> nT (6) 


with |n|?= 1. Then the electromagnetic coupling of fermions which are 
isOspinors reduce to 


AP ,O; (1 + 75) be fi (7) 
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where O, and f; are appropriate y-matrices and functions of A, respectively. 
The two cases (1 + 73) belong to the cases where7 = + 1. Itis to be noted 
that couplings between fields of different » are forbidden so that both ¥, 
and yg should have the same 7. In the d’Espagnat-Prentki scheme [d’Espag- 
nat and Prentki (1956 and 1957)], fields with 7 = +1 correspond to iso- 
spinors of the first kind and those with 7 = —1 to isospinors of the second 
kind. When %, = sg, the transformation has of course to be applied simul- 


taneously to both % and ¥ so that the form is not quite unique in that both 


AOwf, and A~O;7,~f; are invariant under these transformations. How- 
ever if we assume that (7) holds unchanged even in this case, we obtain the 


electromagnetic interaction of an isospinor with itself as A~O; (1 + 74) fj. 
For a spinor field which is an isopseudovector, the analogous transformation 
would be 


E—+ Ef = 83€ 
where 
1 0 0 ¢ 
0, = 0 0 0 E= g° (8) 
0 0-1 = 


with |n|?=1. The coupling A¢O;0,¢f; is evidently invariant under (8). 
This form is also unique provided we assume the cylindrical symmetry of the 
electromagnetic interactions about the third axis in the isotopic spin space. 
This scheme is sufficient to ensure the minimal electromagnetic interaction 
of the baryons apart from the case where a particle is coupled to itself through 
a magnetic moment interaction. A coupling between 2* and p or &- and 
=- can be forbidden by choosing 7 = + i for Z even if we do not assume 
any rotational symmetry for electromagnetic interactions about the third 
axis in the isotopic spin space. 


5. A MOopDEL OF STRONG INTERACTIONS 


Our considerations regarding the mass gauge transformations have 
suggested a model of strong interactions according to which the K-meson 
interactions exhibit ‘cosmic symmetry’ in the sense of Sakurai [Sakurai 
(1959)] and the 2 — 4 —7 coupling is zero, The interaction Lagrangian 
in this model would be 


List. ai lL, + lx 
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ee i —_ - - a er 
= iG \NiystN,7 -L iG, 2 = 2 57 at iGgNay5TN 47 
+ F [(N,N.) K® + (N,N,) K+ + (N,N.) K* 


— (N,N) K® +/.c] (9) 


where, as usual 








7° @° 
N; = (~ ) N,= s- (10) 
with 
Ao — x0 Ae +50 
Y= Sy z= 5; 


We also define 


a+ _ ye | 
N,’ = ee: Ns -( > ). (11) 


which are obtained by replacing A® by — A® in N, and Ns. 


In writing (9), we have assumed equal parity for all the baryons while 
the appropriate matrix | or iy; has yet to be inserted in the K-baryon 
interactions depending upon the parity of the K-meson. 


The Lagrangian (9) involves a totality of four different coupling con- 
stants and does not lead to any result independent of the perturbation theory 
in contradiction with experiments. In order to obtain a further reduction 
in the number of independent coupling constants, we may draw the follow- 
ing analogy between the K- and z-meson interactions and _ the z- and electro- 
magnetic interactions. The 7-meson interactions connect different charge 
states of a field and are known to be symmetric between the different charge 
states. The electromagnetic interactions break this symmetry. This break- 
ing is not arbitrary, the couplings introduced for a triplet for instance being 
such that each member is coupled to itself with a coupling constant equal to 
elz. The electromagnetic interaction of an isotopic triplet thus transforms 
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as the third component of a vector in the isotopic spin space. Since 
in our model, the K-meson interactions which connect different hypercharge 
(U) states of the baryons are supposed to be symmetric and the 7-meson 
interactions are supposed to break this symmetry, it is natural introduce the 
m-meson interactions in analogy with the electromagnetic interactions by 
setting G, = — G,; G,=0 since N, has U= 1, N, has U=—1 and 
N, and Ns have U = 0. This implies that the z-interactions are cylindrically 
symmetrical about the third axis of the three-dimensional hypercharge space 
introduced by Feinberg and Gursey [Feinberg and Gursey (1959)]. The 
interaction Lagrangian now reduces to 


as. ™ iG [Niys7Ni os Nyys7N,] 7 
+ F [(N,N.) K® + (N,N,) K* + (N,N,) K* 


— (N,Nz) K° + h.c]. (12) 


However this interaction is too highly symmetrical and as Pais has 
shown [Pais (1958)] leads to a number of results in contradiction with experi- 
ment. We may break this symmetry by insisting only on the equality of the 
absolute values of the K-meson coupling constants and setting for instance 


faux = fenx = fuax = — Snox = F. (13) 
We then obtain, using (11), 


Lig. = 1G [Ni 7y5N, = N,tysNq]-7 
+ F [(N,Ns) K° + (N,Ns) K+ + (N,N,’) K: 


— (N,N,’) K® + hc). (14) 


Considerations analogous to that of Pais [Pais (1958)] show that this removes 
the 2° — A® degeneracy. However the invariance of (14) under 


N,—>N,, Nz—>Ni, 7—>— 7, A*>— Ae, Kt+—Ke Ke.,.£ (15) 


implies that the nucleon and the cascade are still degenerate. Invariance 
under (15) also implies certain equalities between the amplitudes involving 
the nucleon and the cascade which at present cannot be verified. It is 
interesting to note that the difficulty of introducing the mass difference 
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between N, and N, here is rather analogous to the difficulty of introducing 
the 2+ — 2 mass difference which as pointed out by Gell-Mann [Gell-Mann 
(1957)] is rather too large to be explained as due to purely electromagnetic 
interactions. 


The work of Pais [Pais (1958)] has already shown that such a small num- 
ber of independent coupling constants as two is insufficient to reproduce all 
the known facts. It is however felt desirable to point out the suggestive 
analogy between the z- and electromagnetic interactions on the one hand 
and the K- and z-interactions on the other which does provide us with a 
definite way of introducing the z- and K-interactions. One can of course 
remove the degeneracy of the nucleon and the cascade by choosing different 
values for G, and G; still setting G,=0 or by choosing fan = fonx 
# faux = — /en,x- Such procedures however are rather arbitrary. 


This model of course does not imply that the 7-2-2 or the 7-2-A inter- 
actions are identically null. We always have virtual processes, where for 
instance the 2 dissociates into a nucleon and a K, the nucleon interacts with 
the 7, and the resultant nucleon reabsorbs the K to become a Z or a A, 
These are thus two-step processes and we may perhaps consequently expect 
the 7-2 and z-A interaction cross-sections to be smaller than the z-nucleon 
cross-sections. 


We are grateful to Professor Abdus Salam for his comments and criti- 
cisms in the light of which the paper was revised. He has pointed out that 
as regards the mass gauge transformation, since we have made more use of 
the violation than the validity of the law of fermion mass conservation, the 
conclusions from them should not be pressed too far. 


We also wish to thank Mr. K. Venkatesan for helpful discussions. One 
of us (N. R. R.) is indebted to the Atomic Energy Commission for the grant 
of a Senior Research Fellowship while the other (A. P. B.) is grateful to the 
University of Madras for the award of a Fellowship. 


6. SUMMARY 


The existence of symmetries which can lead to minimal electromagnetic 
interactions have been investigated and three specific invariances to obtain 
this minimal character have been proposed. Applications of these principles 
to strong and weak interactions are seen to lead to certain suggestive results. 
A model of strong interactions arising out of these considerations is also dis- 
cussed. 
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APPENDIX 


The analogy between the 7,- and y;- transformations discussed in Sec- 
tion 4 leads to an isotopic spin scheme for leptons in which y~ and e~ are one 
of the components of doublets whose other components are missing and v 
is a singlet. This scheme is discussed below and it is found that it is able to 
explain many of the features of leptonic interactions in a qualitative way. 


THE Isotopic SPIN SCHEME FOR LEPTONS 


In Section 4, we have discussed the fact that the electromagnetic inter- 
actions of baryons are invariant under the transformation 4—»7y7, with 
|n|?= 1. This is analogous to the y,- transformation of the neutrino field 
and eliminates one of the components of the baryon field in its interaction 
with the electromagnetic field. We now extend the above transformation 
to » and e fields and assume that they possess the following invariance similar 
to the y,-invariance of the neutrino field: 


= — tee 
e = — Tse (Al) 


where the » and e fields are taken to form isodoublets. It follows that yu 
and e are essentially one component isotopic spinors of the form 


‘ oe mee | (A2) 


We further assume that the neutrino is an isosinglet. 


STRONG INTERACTIONS 


It will now be shown that with the isotopic spin assignments given above, 
it is possible to forbid the strong interactions of leptons. We assume 


(a) Strong interactions conserve the total isotopic spin. It is immediately 
clear that w~ and e~ can have no strong interactions since no isotopic scalar 
can be formed involving » or e because of their missing components. 


If we further assume 


(b) Strong interactions conserve parity, v also can have no strong inter- 
actions since it always occurs in the combination 4(1 + yz) »v. 
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It follows that leptons can have no strong interactions. 


In electromagnetic interactions, the combination jf-y- or e-e~ occurs 
which transforms like a scalar plus the third component of a vector which is 
the correct transformation of these interactions for isospinors in the isotopic 
spin space. 


It is interesting to note that only the leptons seem to possess invariances 
under transformations which eliminate one of their components, either in 
Lorentz space (y;-invariance of the neutrino) or in isotopic spin space 
(— 7,-invariance of » and e) while for baryons, such invariances seem to be 
a property only of some of their interactions and not of the fields themselves. 


WEAK INTERACTIONS 


Present evidence in weak interactions seems to rule out the existence of 
neutral lepton currents [Lee and Yang (1960)]. We therefore assume that 


the only leptonic currents relevant for weak interactions are (jv) and (ev). 
These transform as one of the components of an isospinor. We now assume 
that weak interactions involving leptons in which baryons are also involved 


> 
obey the rule |41| = 1/2 or 3/2 while weak interactions involving only 


leptons obey |41| = 0 or |. Strangeness-changing baryonic currents satisfying 
|4S| = 1 necessarily transform as isospinors corresponding to I = 1/2 or 3/2. 


Consequently an interaction of the form J,SNc* J, with |4S| = 1 will 


not obey the selection rule |A1| = 1/2 or 3/2 and are consequently forbidden. 
Thus we get the result that strangeness-changing leptonic decays with 
|4S|= 1 are forbidden. Experimentally, apart from the K—>p + » decay, 
all such decays have low probabilities. The K—»p + v decay is hard to 
explain. Thus apart from this decay mode, we are able to explain at least 
qualitatively the low rates of such processes. 


In this scheme, if |4S| = 2 currents are allowed, we may have decays 
like 3-—+>n +e- +. The K,° — K,° mass difference however seems to 
suggest against the existence of such currents [Okun and Pontecorvo 
(1958)]. 
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SUMMARY 


A quartic profile in terms of the normal distance from the wall 
has been taken and coefficients are evaluated by satisfying one more 
boundary condition on the wall than the usual one. By doing so, the 
limitations about the Reynolds number of the quartic profile adopted 
by Lew (1949) has been removed. The K4rman (1921) Momentum Integral 
Equation has been used to evaluate the various characteristics of the flow. 
A comparative study of Lew’s quartic profile and exponential profile 
together with the quartic profile of the present paper has been under- 
taken and the graphs for the various characteristics of the flow for a 
number of Mach numbers and suction coefficients have been drawn. 
At the end, certain conclusions of general nature about the velocity pro- 
files have been recorded. 


1. INTRODUCTION 


Lew (1949) solved the problem of boundary layer control for Prandtl number 
unity with specific heat Cp at constant pressure as constant and viscosity 
temperature relation given by 4»coT™. He took a quartic profile as an 
approximation to the actual profile in the boundary layer and found that 
his results are applicable up to the Reynolds number 4 x 10° only. To 
overcome this difficulty, he used an exponential profile, which apart from 
other unsatisfactory features, gives substantially higher values of the charac- 
teristics of the flow throughout the entire range of Reynolds numbers. 


In what follows, an attempt has been made to improve his solution 
so that the results are valid throughout the entire range of Reynolds numbers. 
This has been made possible by taking a quartic profile in terms of the dis- 
tances normal to the wall and evaluating the coefficients by satisfying one 
more boundary condition on the wall than the usual ones. 


A comparative study of the present quartic profile with Lew’s quartic 
profile along with his exponential profile has been made and the graphs for 
12 
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the various characteristics of the flow for a number of Mach numbers and 
suction coefficients have been drawn. At the end, certain conclusions of 
general nature about the suitability of the various profiles have been recorded. 


Assumptions.—The analysis is based on the following assumptions: 


: d Pi) 
(i) Ix SF’ 


(ii) Prandtl number oc is unity, 


(iii) (=), = 0, where the suffix ‘0’ denotes the value on the 
plate. 

(iv) Cp = constant 

(v) pooT™ 

(vi) pu = constant. 


Assumption (i) is the usual boundary layer approximation. Prandtl 
number o = (ucp/k) = 1 is taken in order to obtain a particular solution 
of the energy equation in the case of no heat transfer. Actual value of o 
for air is 0-74 and, therefore, it slightly overestimates the viscosity. The 
assumption Cp = constant holds for a wide range of temperature. As for 
the assumption (v), 7 varies from 0-768 to 0-5 depending on the temperature 
range. At room temperature, n is approximately 0-768 but it approaches 
to 0-5 at higher temperatures. In this paper, all the numerical work has 
been carried out for n = 0-768. py = constant is taken only to make the 
velocity profile independent of temperature distribution. 


1. EQUATIONS OF MOTION 


The equations determining a steady compressible laminar boundary 
layer flow on the plate without pressure gradient are the equation of momen- 
tum 


pus + v= 5 (HS), (1.1) 
the meee of — 
2 (pu) + = >, (ov) = 0, (1.2) 
the equation of energy 


) ) _ hee du\? 
pus (cpl) + 90 5 (eT) = 5 (kK) + a(S ’ (1.3) 
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and the equation of state 
p = RpT. (1.4) 
In the above equations x, y are the distances measured along and per- 


pendicular to the length of the plate and u, v are the velocity components 
in the corresponding directions. Other symbols have their usual meanings. 


These equations are to be solved under the following boundary con- 
ditions: 


(i) At y=0, u=0, v = vg = constant (v_, < 0 for suction 
and v, > 0 for injection), 
and (1-5) 
(ii) At y= 8 (x), u=U, T=T,, 
where 6 (x) is the boundary layer thickness at any section x and T, is the 
temperature in the main flow. 


Multiplying (1.2) by uw and adding to (1.3), we get 


dTX oT* i 1) oT 0 = 
wo tee Ges) at a) 6 
where 
Die 
Tx = s YY — u 
“Cp 
when 
e=f, 
T* = constant ..7) 


is thus always a particular solution of the equation (1.6), regardless of the 
distribution of pressure at the free-stream boundary. Since 


dT udu 


oy Cp OV 


vanishes at the wall, the integral (1.7) solves the problem of no heat 
transfer for any given outer velocity distribution. We note that it implies 


Daa 
T+ om u* = T, = constant. 








Tp ed 


Ry teh OE IO tt 
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This particular integral (1.7) makes the analysis much simpler as once 
the velocity distribution is known, the corresponding temperature distri- 
bution can be easily calculated from (1.7). 


The differential equations, governing the variation of u, v, p and T are 
(1.1), (1.2) and (1.6) together with the perfect gas relation (1.4). These 
equations are to be solved under the boundary conditions (1.5). Since, 
the pressure in the boundary layer is supposed to be known and T is given 
by (1.7) in case u is known, therefore, p is known from (1.4). As such the 
solution of these equations reduce to finding the velocity components u and 
v from equations (1.1) and (1.2) with the usual boundary conditions (1.5). 


2. SOLUTION OF THE EQUATIONS 


As a first step towards the solution of the equations (1.1) and (1.2), 
we transform the co-ordinate system from (x, y) to (x, t) plane by the following 
Dorodnitzyn (1942) transformation: 


xX =X, 


y 
— -_"* 
f = few. (2.1) 


The transformed momentum equation in the new variables with the 
assumption (vi) becomes [Kuerti (1951)] 


du _ dU d2u 
Us, +8 y= Ye See em 
where 
” dt 
B=us +», (2.3) 
@ 


and the equation of continuity (1.2) is 


wu , Ow 


The boundary conditions (1.5) take the form 
(i) Att=0, u=0, v = vy = constant 


(ii) Att =&, u=U (2.5) 
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where } 


8 (x) = f “di 


3; 
= fii m inT,) dt. (2.6) 


It is to be noticed that equations (2.2) and (2.4) are of the same form 
as the corresponding equations in the incompressible fluid, the effect of the 
compressibility is confined only in the transformation of the normal co- 
ordinates and normal component of the velocity in the boundary layer. 


Variables are made dimensionless with corresponding quantities at 
the wall in the following manner: 





< ul = ~ 
ee a 2 
i a? p ap Pp 
‘i i ' (2.7) 
= T . be - t = X 
r= = p= , i — ' 
T, la) L L / 
where L is some characteristic length. 
Equations (2.2) and (2.4) in the dimensionless form become 
- vu _ du 1 07% 
" — — — Be” =e? ia 
_—" Re d7? (2.8) 
nr ny) 
se + 0, (2.8 a) 
where 
Ry = lv TeTe 
Ly i 
and the boundary conditions (2.5) reduce to 
A y 
i =% = 0, v = Uo 
where 
_ as <- v q 
t =5¢, “= U, Va. = = 2.9 F 
0 V 2pTy ( ) [ 
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Integrating the boundary layer equation (2.8) across the boundary 
layer from 0 to 8 and using (2.8a) to eliminate 0, we get the well-known 
Von Karman momentum integral equation as: 


5, 
P) — —_" eS vu 
af a (O— a) di — R,o,f = (=), (2.10) 
where 
sS= Rex. 


Since, in the boundary layer problem, we are concerned with the cal- 
culation of skin-friction, which mainly depends on the conditions near the 
body, we should place more emphasis on the boundary conditions near the 
body rather than at the outer edge of the boundary layer. We obtain the 
requisite number of boundary conditions at the body by differentiating the 
boundary layer equation and continuity equation [Howarth (1934)]. By 
doing so, the difficulty experienced by Lew in the case of his quartic profile 
is eliminated. These conditions have been first used by Torda (1951) for 
the incompressible case. 


We assume the following form of the velocity profile: 


= a(s)7t + b(s) 7? + c(s) tr? + d(s) A, (2.11) 


CHS 


where 





a, (#) = 1 (2% 
at 7=0 0 dt 0 Re m),° 
= (5) — | (2% 
. or? . Re 7), ° (2.12) 
and 
a =U, 
at r= l 7} 
<= 0, 
or / 


A2 














| 
| 
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give the following values of the coefficients: 











ia 24 es 12M 
(= mapMais: OmepeMaie | 
ion d() = 3+ M+) O19) 
=P eM +18’ °- MRF OM +18” 


where 
M = Redo dt. 


It is clear that for no real value of M, the coefficients become infinite. 


Substituting the expression for the velocity profile (2.11) with the co- 
efficients (2.13) into the Von Karman Integral equation (2.10), we get 


Ren,” 
o dy 





4 M (MS + 18M + 157M4 + 914M® + 3588M? + 7572M + 7308) 
35 (M2 + 6M + 18)?(M3 + 6M? + I8M + 24) 





e) 


(2.14) 


Integrating (2.14) between the limits 0 to M and evaluating the constant 
of integration from the condition: 


M=0 at s=0. 


we get 


on (R) s*) = - {1 2 eet 1 M..*}"" 


- -4M2-+.99M 
= {M T (M2 + 6M + 18) 





2 
— 3log M+ OM + 1S _ 5.0848 log tm 





2 ° . 
+ ong 





M 
M +6 





+ 17 tan- 





2- 5088 
oa 1 
17.4 tan een O13} ° (2.15) 


rah tae ig oc an 


Teen a RI AES A TN 


erin Nt. OE APN CO IR 





wht 


an 
pre 
the 
(2. 














Compressible Boundary Layer on Flat Plate 19 


where 


M.. = 2 a being the velocity of sound in the main flow. 

The equation (2.15) establishes the relation between the parameter M 
and the actual length of the plate €. The frictional drag coefficient, velocity 
profiles and boundary layer thickness can now be calculated in terms of 
the parameter M and related to the actual run of length é through equation 
(2.15). 


The coefficient of average friction is defined by 


s ( 4) - 














Cr = ¥p,U'x (2.16) 
— nr 
=-— { pet .} F(M) 
Vo 2 
(Rex) (79 
“ 2(7 , (2.17) 
a * Meo? 
where 
F(M) = & , M(Mt + 13M® + 89M? + 278M + 406) 
— 35° (M* + 6M + 18)? 
and 
U - 
Ree = >> 


The relation between the transformed ordinate ¢ and the physical ordinate 
y can be established with the help of the relation (2.1), viz., 


y=EL f (1 — adr. (2.18) 
0 
Using (2.11) with (2.13) in (2.18), we get 
Voy _ n+1 
ny ~'} 


y—! 2 _G(M,7) 
x ets 1. (2.19) 








1+7>M, 
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where 


G(M, 7) = {192 + 144Mr + + M272 


r3 
(M2 + 6M + 18)2 
4 8r3 [2M — 3(M2-+4 2M + 2)] 
+ 8/7 Mr4 [2M3 — 9(M2 + 2M + 2)] 
— 3 M2r5 (M2? + 2M + 2) 


+ 7 (M? + 2M + 2)*}. 


Putting t= 1 in (2.19), we get the following equation giving us the 
boundary layer thickness in terms of the parameter M: 


Vy 


xM {1475 Mr — GM, 1}, (2.20) 


where 
2 | 5 
35 (M2 + 6M + 13)2 ¢ 


+ 1460 M + 2590). 





G (M, 1) = M*‘ + 65M + 431M? 


Lastly, the velocity profile, given by (2.11) and the temperature 
distribution, given by (1.7) are related through equation (2.15) and (2.19) 
for any Mach number M,, and suction coefficient v9/U. 

3. For sake of comparison, we record below the € — A relation, C,, 
y—t felation and thickness of the boundary layer as obtained by Lew* 
where his A is equal to 4M of the present paper 


9% y—1 mr 30 | 
€ =(1+7> Mo) [M@-yyatle 
9, (+3) Bo. is 
+s"srt uss ye FAFDH “ 


* The expressions (3.1), (3.2), (3.4), (3.7) and (3.8) in Lew’s paper contain some misprints. 
For example, the log and arc tan terms in (3.1) are wrongly printed as 


(a + 3% 


08 309+ AS)” 


arc tan {A ¥3(A + 2)}. 











the 
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Cy = =m 2)" f(A) 





—2(% (3.2) 
id Y— Moo 2 











where 
x A + 26A + 37) 
fO) = 3 ee 
Also, 
a Voy _ nN+1 
Vy — ) 
— | G (A, 7) 
x A wiaUcouet M_2 
— Feel arp (3.3) 
where 
G (A, 7) = WS $3 [840 + 1260Ar + 168 (3A? — 8A — 6) tr? 
— 140 (8A? + 3A — 3) r? + 40 (25 A2 + 33A + 9) 
— 105 (4A? + 7A + 3) + 70(A 4 1)? 7°], 
and 
= Vo (xX) _ y= 1 2)” 
—— sill 5 te 
x {1+75—MQt fl — g QI}, (3.4) 
where 


__ 3402 + 221A + 367 
8) =a + 3 





Lew’s velocity profile gives values of u/Us> 1 for values of A< — 1°5 
so much so that when A = — 3 the velocity becomes infinite. The graphs 
for U8/v, against € shows abrupt changes when A < — 1-5 showing thereby 
that his expression for boundary layer thickness does not hold for values of 
A< — 1-5, We can see that this corresponds to Reynolds number Ux/», 
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= 4x 105. Hence, the above results are applicable only for Reynolds 
number < 4 x 10°. 


Lew overcame this difficulty by choosing the following exponential 
form of the velocity profile 


5 =1—e7[l — rk(s)], (3.5) 


where k(s) is so chosen as to satisfy boundary condition 
=6(%),= G2) 
Vat), ReNdi? Jo 

on the plate, so that 


1+ 2A 


k(s) = — roe (3.6) 


NI — 


It is clear from (3.8) that this velocity profile satisfies the boundary 
conditions 


ws, dw ae ot _ 6 
oo ae a | ee 
at r= oo and not at r= 1. 


For exponential velocity profile, Lew finds that 


=(B) ()=5f1+ 24 " 
+ ()'()= [eM 








A (3A-+5) l A 
and 
-f ae —y— a, gl™  A(I2A2 + 24d + 11) 
Cy mmnages f+ a Mt cae 
(Rex) (73 
. fs : (3.8) 


4. DISCUSSION OF THE RESULTS 


In Fig. |, the velocity profiles for M,, = 1 and 3 have been drawn. It 
is shown that the effect of Mach number on the velocity profile is similar 
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Fic. 1. Velocity profiles variation along plate length. 
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Fig. 2, Boundary layer thickness variation with suction coefficient 
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Fic. 3. Boundary layer thickness variation with Machnumber along plate length. 
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to the effect of a favourable pressure gradient along the length of the plate. 
With the increase of Mach number, the profiles become nearly linear 
[Karman and Tsien (1948)]. In Fig. 2, graphs for the boundary layer thick- 
ness against the length of the plate € for M,, = | and the suction coefficients 
v»/U = — -001, — -002 and — - 003 have been drawn. It is found that 
boundary layer thickness decreases markedly as the suction coefficient in- 
creases and the present profile gives lower values of the boundary layer 
thickness than the Lew’s quartic profile. In Fig. 3, it is shown that as 
Mach number increases, the boundary layer thickness increases. In Fig. 4, 
local shearing stresses as given by the quartic profile of Lew and of the 
present paper have been compared and it is found that the present profile 
gives lesser values of the local shearing stresses for M,, = | and (v»/U) = 
— -001 for Reynolds number up to 8-5 x 10*. Beyond this range, the two 
curves are the same in the range of validity of Lew’s profile. In Fig. 5, 
various curves for the drag coefficient of the present profile for M,, = 1 
and various values of the suction coefficients have been drawn and it is 


O10 
oo9 





3 45 67869 2 3 45 678 2 3°45 67869 - 3.64865 


Ux 
Rx s \- ae 
Fic. 5. Variation of friction coefficient with suction. 


found that, as the suction coefficient increases the drag coefficient also 
increases. This is expected as due to suction the boundary layer thickness 
decreases markedly as shown in Fig. 2, therefore, the curve for the velocity 
profile will rise steeply to gain the value of the velocity in the outer flow. 
As such, 7 =p,(du/dy), will increase and therefore, the drag co- 
efficient will increase. However, due to the thinning of the boundary layer 
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there will be increased stability and the velocity profiles will be less prone 
to transition. It is, therefore, expected that the net drag of the body will 
be reduced significantly due mainly to delayed transition. 


In Fig. 6, it is shown that as the Mach number increases, the drag 
coefficient decreases. Also, a comparative study of the drag coefficients of 
the present profile with Lew’s quartic profile and exponential profile for 
M,, = 0, 1 and (v,/U) = —-001 has been made. It is found that in the 
range of applicability of Lew’s quartic profile, the present profile gives slightly 
lower values of the drag coefficient up to a certain Reynolds number and 
then, the two curves coincide. For M,,=1 and v,)/U = — -001, the 
exponential profile gives sufficiently higher values of the drag coefficient. 


Comparing the results of the three profiles for the particular case 
M,,, = 9 and (v,/U = 0) with the Blasius profile, we find that Lew’s exponen- 
tial profile overestimates the local shearing stresses by 19-2 percent. and 
his quartic profile overestimates by 3-31 per cent. while the present quartic 
profile underestimates the local shearing stresses by 6°91 per cent. 
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Fic. 6. Comparison of skin friction coefficients. 


5. CONCLUSIONS 


A comparative study of the three profiles have been made. The present 
quartic profile satisfies two boundary conditions at the plate and two at 
the top of the boundary layer, Lew’s quartic profile satisfies one boundary 
condition at the plate and three at the outer edge of the boundary layer and 
the exponential profile satisfies one boundary condition at the plate and 
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infinite boundary conditions at + = (t/5s)= oo. The following conclu- 
sions have been derived: 


(i) Lew’s quartic profile hold in the limited range of Reynolds number 
while the results of our profile and exponential profile hold in the entire 
range of Reynolds number. 


(ii) Characteristics of the flow, given by the exponential profile are 
considerably higher than the present quartic profile or Lew’s quartic profile 
while these characteristics, given by our profile have lower values as com- 
pared to Lew’s quartic profile up to a certain Reynolds number in the range 
of its applicability. 

(iii) Nature of the characteristics of the flow is the same in all cases 
as we get in the case of impermeable wall, namely, with the increase of Mach 
number, the boundary layer thickness increases, the drag coefficient decreases 
and the velocity profiles become almost linear. 


(iv) As the suction coefficient increases, boundary layer thickness 
decreases. This thinning of boundary layer results in delayed transition. 
Therefore, although the drag coefficient increases as the suction coefficient 
increases, the net drag on the body will be reduced due to delayed transi- 
tion. 
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INTRODUCTION 


IN earlier parts of this series were reported the results of studies on a number 
of aliphatic esters and in particular of ethyl acetate in a variety of solvent 
compositions. Specific solvent influences were noticeable in which dielec- 
tric constant was only one of the factors and abnormal features were noticed 
where the organic solvent component in the aqueous system was a small 
fraction of the mixture. It appeared necessary to study some aromatic 
esters as well as the esters of polybasic acids. The present report presents 
the results of studies with some benzoic esters. 


EXPERIMENTAL 


The solvents used were purified as in earlier parts while the esters were 
invariably freshly purified and distilled before use. The physical con- 
stants of the esters are given below: 


Ethyl benzoate b.p. 212-9° C. Np” 1-5028 Dipole moment 1-95 
o-Nitro Ethyl benzoate b.p. 173-0/18 mm. ip 1-5211 ib " 3-79 
p-Nitro Ethyl benzoate m.p. 57° i 3-84 
Benzyl benzoate b.p. 323° Np” 1-5641 - ‘i 2-08 


(All the physical constants were measured using the samples purified 
for the experiment. Dipole moment values refer to benzene solutions.) 


The course of the reaction was followed by standard procedures, ana- 
lysing aliquot parts. Freshly standardised dilute ammonia solutions were 
used for estimating excess acid used for chilling the reaction mixture and 
a screened indicator composed of neutral red and methylene blue was used. 
Temperature control was accurate to 0-02° C. in the range studied and was 
periodically checked. 
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DISCUSSION OF RESULTS 


Unlike the studies with ethyl acetate, the range of solvent compositions 
that could be used was limited by solubility of the reactants. In every 
instance, strict second order rate constants could be got and the graph of 
log (a — x)/(6 — x) against time was linear within the limits of experimental 
error (Fig. 1). The relevant rate constant data are presented in Table I. 
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It will be noticed that in all the instances, the reaction rate decreases 
with a decrease in dielectric constant, indicating repulsive interactions of 
the ion-dipole system if one is to assume what is not true, that free energies 
of hydration of the ion as well as of the transition state are unimportant. 
At the concentrations of reactants used, the Debye-Huckel limiting law is 
applicable and no serious error is introduced in applying the Kirkwood 
approximation.‘ 

In dealing with mixed solvents, one has to choose between iso-com- 
position and iso-dielectric conditions for comparisons. Svirbeley and 
Warner! object to the use of isocomposition data as such and use the relation. 

log dD 
Ey — Eapp + R| Sp |. ry 
»(z) 
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TABLE I B 


The alkaline hydrolysis of benzyl benzoate in acetone-water system 








Temperature ° C. 











Solvent 30 35 40 
composition 
Acetone% v/v Conc. of k, Conc. of k, Conc. of ks 
reactants x 108 reactants x 108 reactants x 108 
x 10° M x 10°M x 108 
50 4-824 9-31 4-926 13-40 5-209 19-04 
60 5-226 7-95 5-014 11-57 5-129 16-50 
70 5-098 5-901 5-024 8-70 4-999 12-69 











Note.—The concentration of the hydroxide used was practically the same as that of the 
ester except those marked * where the OH’ concentration was about 0-001 M higher. All rate 
constants are in units of moles/litre, sec.-* 


The degree of linearity in the graph of the reciprocal of temperature against 
the logarithm of the rate constants in the present studies indicate that appa- 
rently the correction term is not appreciable and we are justified in using 
the more convenient isocomposition data in Table II. 


Laidler and Eyring® have shown that in the alkaline hydrolysis of ethyl 
benzoate in alcohol-water systems, the empirical relationship reported by 
Evans and Polanyi? TAS=a4H+8 holds good. While our results 
show clearly that the trends in both 4H and JS appear to be similar, the 
expected linear relationship is not satisfactory. Apparently, where non- 
hydroxylic solverts are present, the solvation of the reactant species as well 
as the transition state are influenced to different degrees by the dielectric 
constant in the microscopic region as well as the topology of the transition 
state. The results, however, bring out the Fairclough-Hinshelwood rela- 
tionship between the two Arrhenius parameters.* Even here, we have to 
treat each compound separately in a given solvent pair and changes in solvent 
pair systems also alter the individual values though the trend for a given 
compound may be maintained. 


For an ion-dipole reaction, assuming a solvent continuum, it has been 
shown that the rates and activation energies involve functions of the dipole 
moments of the reacting species. The reaction centre in these esters is the 
COOR group and substituents in the benzene ring can modify the group 
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TABLE II 
Ester Solvent E 4H logy, 4S 
dioxan vol.% __ k.cals. k.cals. PZ e.u. 
Ethyl benzoate oi 40 12-220 11-608 6-800 —27°51 
50 12-760 12-148 7-300 —25-17 
60 13-230 +=12-618 7-600 —23-84 
70 13-950 13-338 7-760 —23-11 
o-Nitro ethyl benzoate 40 10-800 10-818 6-696 —27-97 
50 11-260 10-648 6-908 —27-01 
60 11-610 10-998 7-070 —26-26 
70 12-010 11-398 7-293 —25+24 
p-Nitro ethyl benzoate 60 11-570 10-948 8-102 —21-54 
70 12-280 11-668 8-552 —19-47 
Benzyl benzoate ie 50 12-470 11-858 7-066 —26:27 
60 12-790 12-178 7220 —25-59 
70 13-280 12-678 7-508 —24-25 
Acetone water 
system % 
Acetone by 
vol. 
50 13-550 12-938 7*737 —23-21 
60 13-960 13-348 7-964 —22-17 
70 14-530 13-918 8-250 —20-91 





moment by altering the apparent charge on the carbon of the carbonyl. 
The introduction of the nitro group can then be expected to reduce the 
activation evergy for the reaction as is actually found to be the case. The 
position in the benzene ring does not appear to make any serious difference 
to the activation energy but as may be expected, there is a pronounced dif- 
The value is considerably higher with 


ference in the entropy of activation. 








a. in -~ my 
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the ortho-nitro compound. The proximity of the negative end of the N-O 
dipole hinders the approach of the OH’ ion from certain directions. Altera- 
tion of the alkoxy group is not without effect and the replacement of the 
ethyl by benzyl leads to a slight reduction in the activation energy following 
the trend of the dipole moment. 


Turning to solvent influence on the reaction, both rates and activation 
energies follow the trend in the dielectric constant with rates decreasing and 
activation energies increasing as the water content decreases. 


It has been observed with ethyl acetate hydrolysis that abnormalities 
are noticed only for small additions of the organic solvent to water while 
in the region of compositions used in the present study, the characteristic 
correlation of Kirkwood‘ between Dielectric constant and rates holds good. 
The one compound studied in both dioxan-water and acetone-water system 
shows again the differences noticed in earlier parts. The reaction in acetone 
is consistently slower. The activation energy as well as the entropy of activa- 
tion are also much higher and beyond the limits of experimental error. The 
difference has to be traced to the specific property of acetone in forming 
a solvated molecule of the type 


OH, ,H,O 


‘ 0° 
Hy AN y 
0---H,C CH,---6 
HZ \u 


so that the system as a whole has to use also the extra energy for rupture 
of the bond between water and acetone in addition to the normal energy 
needed for the reaction. The increase in value of the entropy of activation 
is also consistent with this structure breaking part of the mechanism. 


SUMMARY 


The alkaline hydrolysis of some benzoic esters in dioxan-water mixtures 
as well as that of benzyl benzoate in acetone-water has been studied. The 
observations indicate a group influence consistent with the ease of formation 
of the transition state by interaction between the hydroxyl ion and the 
carbonyl group of the ester in the rate determining step. The influence 
of dielectric constant is what might be expected for an ion-dipole reaction 
but specific solvent influences cannot be excluded. 
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IN a paper Brij Mohan® has given the following theorem connecting different 


classes of self-reciprocal functions: 


If f(x) is R, and 


k+ico 
P(x) = 55, f APE + de + PE + by + 45) 060) 8d, 
ate (1.1) 
where 
O<k<1, and w(s)= (1 —5), (1.2) 
then the function 
g(x) = i P (xy) f(y) dy (1.3) 


is R,. 
Here we proceed to take up the converse problem, namely ; 


If f(x) is R, and g(x) is R,, what will be the character of the function 
P(x) in (1.3)? 


The solution of this problem is covered by the following: 


THEOREM 1: In (1.3) if f(x) is R, and g(x) is R,, then P(x) will be 
given by (1.1). 


Multiplying (1.3) by x? and integrating with respect to x between 
the limits 0 and oo, we get 


T xg (x)dx = f x®ldx f P(xy) f(y) dy 


ai 


f(y) dy J x8-1P (xy) dx, 
35 
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provided the change in the order of integration is justified. Therefore, ye 
have 


0 
Since f(x) is Ry, we have? 


6+ico 


S() = 55, f HLA + du + 45) 4 (9) xs, 


where 0< c< 1 and ¢#(s) satisfies (1.2). 
Hence 

Ff (e) x tde = T(E + du + 49) 969). 
changing s into 1 —s5, 

fx8f(x)dx = 2-7 Q + te — 45) (5), 
since 

#(s) = d%(1 — 5). 


Similarly, as g(x) is R,, we shall have® 
f x8-lg (x) dx = 2*P'(4 + 4v + 45)X (s), 
0 


where X (s) satisfies (1.2). 


Hence (1.4), with the help of the above known integrals, is easily reducible 
to the following integral 


co 


| 28 T'($ + 4 + 45)X (s 
f us ip (u) du = ral are * > 
= 2P(4 + 44 + 45) (4 + 40 + 4s) gd, (9), 
where 


— X (s) 


py (Ss) = arg a 





thy +45) TE 4+ hu —I)F) 


S g(x) x*"dx me vy S f(y) dv § uP (u) du. (1.4 F 
0 ° 
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and hence satisfies the relation 
b, (s) = 4, (1 — 5). 


Therefore, by Mellin’s Inversion formula,! we have 


k+ico 


Pibe ns f 2°P(4 + dy + 45) P(4 + dv + 45) 4, (8) w-8ds, 


k—ico 
where #, (s) satisfies (1.2), which is the same as the function given by (1.1). 


The generalisation of the above theorem, by the same author’ is given 
by the following: 


If f(x) is of the form 


e+ico 


f(x) = = f 2*#P(4 + $u + 45) P(¥ + 4v + 45) (3 + 4A — 4s) 
x o(s) x-Sds, (1.5) 


where #(s) satisfies (1.2), and P(x) is given by (1.1), then the function 
g(x) given by (1.3) belongs to R,. 
The corresponding converse result will be: 


THEOREM | (a): In (1.3) if f(x) is given by (1.5) and g(x) is R,, then 
P(x) will be given by (1.1). 


This can be established in the same way as Theorem 1. Another prob- 
lem of a similar nature is: 


If g (x) is R, and P(x) is of the form (1.1), what will be the character 
of the function f(x)? 


THEOREM 2: In (1.3) if g(x) is R, and P(x) is given by (1.1), then 
the function f(x) will be of the form 
e¢+ico 


f@=x, f 2 (4 + du + 45) (5) x-8ds, 


where %(s) satisfies the relation (1.2), and hence will be R,. 


The corresponding converse for the generalised theorem referred tg 
above would be; 
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THEOREM 2(a): In (1.3) if g(x) is Ry and P(x) is of the form (1,1), 
then f(x) will be given by (1.5). 


The above theorems could be established in the same way as Theorem |, 


2. Another theorem due to the same author® states: 


If f(x) is R, and 
k+ico 


P(x) = 50, f TOE + du + ATR + dy — b5) (9) 8s, (2.1) 


k—ico 


where 0< k < Il, and w(s) satisfies (1.2), then the function 


g(x) = f P(y) f(xy) dy (2.2) 


belongs to R,,. 


The converse theorem would be: : 

THEOREM 3: In (2.2) if f(x) is R, and g(x) is R,, then P(x) will be 
given by (2.1). 

In Theorem 3 if we put » = », it will give rise to the following simple 
Corollary : 

Corollary.—in (2.2) if f(x) is R,, and g(x) is R,, then P(x) will be 
given by 

k+ico 


P(x = f A(s) x-*as, 


k—ico 
where ‘A (s) satisfies (1.2). 


In another paper Brij Mohan’ generalised the above theorem in the 
following form: 


If f(x) is of the form 


e+400 


f) = 5, f BPE + e+ dT] + be — 45) 


@e-ico 


x F'(4 + 4A + 45) X (5) x-Sds, (2.3) 


where X (s) satisfies (1.2), then the function given by (2.2) is Ry, provided 
that P(x) is given by (2.1). 
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The corresponding converse theorem will be: 


THEOREM 3(a): In (2.2) if f(x) is of the form (2.3) and g(x) is R,, 
then P(x) is given by (2.1). 


When p = »v, in (2.3), the corollary given under Theorem 3 also holds. 
We have two more similar theorems: 


THEOREM 4: In (2.2) if g(x) is R, and P(x) is of the form (2.1), then 
f(x) will be given by 


o+ico 


SO) = a5; f MIE + te + 49) 460) x8 ds, 


o—ico 


where %(s) = %(1 —s), and hence will be R,. 
If 1 = v, in Theorem 4, we get the following Corollary: 


Corollary.—In (2.2) if g(x) is R, and 


k+ico 


P(x) = ms f A(s) x-8ds, 


k—i{Cco 


where A(s) = A(1 — 5), then the function f(x) will also be R,. 


THEOREM 4(a): In (2.2) if g(x) is Ry and P(x) is of the form (2.1), 
then f(x) will be given by (2.3). 


When p» = »v, a corollary similar to that given above will be obtained. 


The proofs of Theorems 3, 3(a), 4 and 4(a) are similar to that of the 
theorem established above. 


3. In another paper the same author® has given the following theorem. 


If f(x) is R,, the function 


g(x) =4 f Q (log) £0) dy 
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is R,, provided that 


k+ico 
Q(x) = <7 f I ($+4y+-4s) P'(3+4v—4s) w (s) e8ds  (x>0) 
k--ico 
wil (x<0), 
(3.2) 
where k > 0 and w(s) satisfies 
w(s) = w(l —s). (3-3) 


The converse problem here would be: 


If f(x) is R, and g(x) is R,, what will be the character of the function 
Q(x) in (3.1)? 


The solution of this problem is covered by the following: 


THEOREM 5: In (3.1) if f(x) is R, and g(x) is R,, then Q(x) will be 
given by (3.2). 


Substituting u for log(x/y) in (3.1), we get 
g(x)= f Q(u) fre u) edu. 


Multiplying this equation by x*-', and integrating with respect to x, bet- 
ween the limits 0 and oo, we have: 


fe (xs) 8 *de = i x3 “dx ,Q (u) f (xe") eMdu. 
° 0 ° 


Let x = ye. 


Substituting this value for x on the R.H.S. and rearranging the terms 
we get, 


f x81g (x) dx = JS f(y) yS"'dy ff Q(u) edu, (3.4) 
0 0 


0 


provided the change in the order of integration is justified. Given that 
the functions f(x) and g(x) are R, and R, respectively, we have® 


co 


J x8"f (x) dx = 2*F' (4 + 4u + 45) 4 (5), 


) 
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and 


f x81 (x) dx = 2#(4 + 4y + 45) X (5), 


where both #(s) and y(s) satisfy the relation (1.2). 


Substituting these values for the respective integrals in (3.4), we get 


~ ans, UG + dv +49)X (5) 
s —— Pr (4 + du + 45) d(s)° 





changing s into (1 — s), 


Fee PS — KG) 
JS werd = FETE HO. 


= (4+ hy + 4s) F(R + dv — $5), (5), 








where 
ee - 
7 P(g + 4p + $5) (3 + Su — 45) o(s)’ 
and hence satisfies the relation #, (s) = 4, (1 — s). 


Hence, by applying another type of Mellin’s Inversion Formula,? we 
obtain 
k+ico 


Qa) = ah, J PE + de + ALE + by — 45) th (6) eFds (x > 0) 


k—ico 
=f (x < 0), 
which is the same as (3.2). 
In another paper* the present writer has generalised the theorem 
due to Brij Mohan in the following form: 
If f(x) is of the form 
| ¢+iCo 4 
$0) =, f BPE + wt 49 PQ + dy — 45) 
6400 


x P(A + 4A + 45) b(s) x-Sds, (3.5) 
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where 0 < c < | and ¢(s) satisfies the relation (1.2), the function g (x) given 
by (3.1) belongs to Ry, provided that (3.2) holds. 


The corresponding converse result for this theorem would be: 


THEOREM 5 (a): In (3.1), if f(x) is of the form (3.5), and g(x) is Ry, 
then Q(x) will be given by (3.2). 


The proof of this theorem is similar to that of the above theorem. 
Another problem of a similar nature is covered by the following. 


THEOREM 6: In (3.1), if g(x) is R, and Q(x) is given by (3.2), then 
f(x) is Ry, 


and the corresponding result for the above generalised theorem would be: 


THEOREM 6(a): In (3.1), if g(x) is Ry and Q(x) is of the form (3.2) 
then f(x) will be given by (3.5). 


The proofs of these two theorems run parallel to that of Theorem 5 
already proved. 


4. Another theorem due to Brij Mohan® is as follows: 


If f(x) is Ry, the function 


g(x) = f Q(los 2) FO) a 4.1) 

is R, provided that 
QW =5, f 28P (4+ 4dut-4s) P(d+4v-+45) w(s) eds (x > 0) 
mit (x < 0), 


(4.2) 
where k >0 and w(s) satisfies the relation (1.2). 


The converse problem is covered by the following 


THEOREM 7: In (4.1), if f(x) is R, (R,) and g (x) is R, (R,) then Q (x) 
is given by (4.2). 


The present writer* has generalised the above theorem of Brij Mohan 
in the following form; 
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If f(x) is a function of the form 


fo) =x, f MPA + +49 TA+ b +49) 
x I'(¢ + 4A — 4s) bb (s) x Sds, (4.3) 


where 0< b< 1, and ¢(s)=y¥(1 — 5), 
then the function given by (4.1) is Ry, provided that Q (x) is given by (4.2). 


The corresponding converse result to this theorem would be: 


THEOREM 7 (a): In (4.1), if f(x) is of the form (4.3), and g(x) is 
Ry, then Q(x) is given by (4.2). 

Two more converse results of a similar nature corresponding to Theorems 
6 and 6 (a) for the above-mentioned theorems of Brij Mohan and the present 
writer would be: 


THEOREM 8: In (4.1), if g(x) is R, and Q(x) is given by (4.2), then 
f(x) is Ry. 

THEOREM 8(a): In (4.1) is g(x) is Ry and Q(x) is of the form (4.2) 
then f(x) will be given by (4.3). 


The theorems of this article can be proved in the same way as those 
of article 3. 


I offer my thanks to Dr. Brij Mohan for his help and guidance in the 
preparation of this note. 
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INTRODUCTION 


Tue T.I.F.R. computer is a parallel, binary, asynchronous machine. It is 
controlled by a stored programme of single address instructions and has a fast 
access ferrite core memory, consisting of 1024 locations. The input to the 
computer is by means of a punched paper tape (5-hole commercial teletype 
tape) and the output can be either printed out directly or punched on a paper 
tape again. 

In an earlier paper,’ the principal features of the engineering and system 
design of this computer were dealt with. In this paper we shall consider 
in some detail the logical and circuit design of its arithmetic unit, control 


and memory. While some knowledge of the logical design and organisation 
of digital computer is assumed on the part of the reader, this paper is self- 
contained as regards the design of the T.I.F.R. machine. 


THE ARITHMETIC LOGIC AND THE CIRCUITRY 


In the T.I.F.R. computer, all the internal operations take place in the 
binary notation and numbers are stored in the main memory as signed 
binary fractions restricted to the range (—1, 1). A number in its normal form 
consists of a sign (which is 0 for + ve numbers and a | for —ve numbers) 
followed by 39 binary digits (bits for short) which give its absolute value as 
a binary fraction. The binary point is understood to be fixed immediately 
after the sign bit. The complement form of a number is obtained by changing 
all its 1’s to 0’s and vice versa, the sign bit is left unaltered in the process. 


A schematic diagram of the arithmetic unit is shown in Fig. 1. It consists 
of a memory register, two (double rank) shift registers (the Accumulator and 
M.Q. register respectively), an adder and a pair of true/compl. gates. The 
contents of the accumulator are presented to the adder either in the true or 
the complement form, depending on whether the operation called for is addi- 


* Now on study leave at Brown University, U.S.A, 
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tion (Multiplication) or subtraction (Division), respectively. In some cases 
the result of a subtraction needs to be recomplemented. To achieve this, 
facilities have been built into the adder so that its output can be gated into 
the accumulator either in its true or complement form. 
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Fic. 1. Schematic Diagram of Arithmetic Unit Asymmetrical. 


The three registers, Memory register, Accumulator and M.Q., are essen- 
tially storage units for holding the numbers involved in the arithmetic opera- 
tions. The basic flip-flop used is of the asymmetrical type as shown in Fig. 2. 
By definition the flip-flop represents the binary digit 1 when the left hand 
triode is cut off and right hand one is conducting. As seen from Fig. 2 the 
output voltage at P for ‘1’ and ‘0’ is + 30 and —20 volts respectively. 
The transfer of information is achieved by first clearing, i.e., making the 
flip-flop ‘0’ and then gating. Clearing is done by dropping the plate supply 
(+ 110 Dr) from + 110 V to about + 40 V. The T.I.F.R. machine being 
parallel, it is necessary for all the 40 flip-flops of a register to be cleared 
simultaneously. This is done by giving the + 110 supply through a driver 
whose output goes down to + 40 V when pulsed. Gating is done by applying 
a + ve pulse rising from —100 V to —10 V. Depending on whether the 
input is ‘ 1’ (+ 30 V) or ‘0’ (—20 V) the flip-flop changes to the state ‘1’ 
or ‘0’ respectively. The clear and gate drivers are discussed in detail later on. 
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As shown in Fig. 1, the accumulator is a double rank (41 bits) register. 
The extra bit is the overflow bit (cf. below). It consists of an Upper Register 


12aT? 6. 42aT 
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Fic. 2. Flip-flop. 


and a Lower Register. Provision is made for shift UP (the contents of the 
Lower Register to the Upper Register), shift Dn (the contents of the Upper 
Register to the lower Register), shift left (the contents of 2-‘"+) Upper 
Register to 2-" Lower Register) and shift right (the contents of 2-‘"-)) Upper 
Register to 2-" Lower Register) and gate adder to the Upper Register. 

M.Q. is also a double rank 41 digit Register capable of shift UP, shift 
left, shift right and gate from Memory Register to Lower Register. The 
extra bit is placed as the 2-®° bit. It serves as the sense bit. 

The adder is of the fastlogical type and incorporates in itself special 
carry by-pass facilities* to speed up the carry propagation time. In an adder 
the factor responsible for serial propagation of carry is given by the 
logical equation 

Cn = Cn ° (An + Bn) 
where A, and By are the inputs (corresponding to 2-" bit) to the adder from 
Memory Register and Accumulator and Cy_, is the carry from the previous 
stage. The above expression means that Cy, can be formed from Cy_, 
if either of the digit inputs is present. Rewriting the above equation 


Cn = Cn-a (An + Bn) 
= Cy. (Ana + Bai) (An + Bn) 


os eer (An + Bn). 


The carry output Cp, from a lower digit stage P, can cause Cy, in case at least 
one dilgit input in all the intermediate stages is present. If means are provided 
for Cp to sense in parallel the digit inputs of all the intermediate stages 
@+h @+H .«....-:. n, Cp will cause Cy without any delay. 
One can achieve this for 7 varying from 2 to N and p from 1 to n — 1 (for an 
N stage adder), in which case no carry delay should occur. This, however, 
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is very cumbersome for N = 40 (the word length in the computer), considering 
the amount of electronic circuitry involved. But a compromise can be reached 
as in Fig. 3 which shows the carry by-pass logic and Fig. 4 which gives the 
circuit diagram as incorporated in the T.I.F.R. computer. 
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Fic. 3. Carry By-pass Logic. 
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Fic. 4. Carry By-pass Circuit. 





A number of 5 adder units are grouped together. There are 8 groups, /.e., 
8 carry by-pass units, the last group having only 4 units. Within this 
group carry propagation is serial, whereas a carry input incident on the least 
significant stage of the group can pass on to the next group directly in case 
at least one of the two digit inputs of every stage in the group is present. The 
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maximum carry delay occurs when carry is generated in the first stage of the 
first group, and dies down in the last stage of the last group, after passing 
through 4 stages, by-passing 6 groups and again going through 3 more stages, 
Assuming that one by-passing introduces the same delay as that of a single 
digit stage, the maximum carry delay is 13 times the carry delay per stage. 
In the present adder the maximum carry delay is 3p sec. 


All proper arithmetic operations should result in numbers which again 
fall within the range (— 1, 1). But there are no built-in checks to detect 
improper arithmetic operations. An extra bit, called the overflow bit, has 
been provided at the left end of the accumulator. In the case of genuine 
overflow (i.e., addition of numbers with the same sign, or subtraction of 
numbers with opposite sign), it is fed into the overflow bit and an overflow 
indicator is set to 1. A “jump on overflow ” order enables the programmer 
to make use of this information in controlling the details of his computations. 


THE CONTROL LOGIC AND THE SEQUENCING UNITS 


As has already been mentioned, the T.I.F.R. computer is of the 
asynchronous, single address type. The details of its control logic are effec- 
tively determined by these two features. A coded instruction in the machine 
requires 20 bits for its specification (the first 9 bits for its order part and 
last 11 bits for its address part). Because the word length in the memory 


is 40 bits, instructions are stored in pairs in each memory location. They 
are identified as the left hand (L.H.I.) and right hand (R.H.I.) instructions 
in that location. In magnetic drum orders two addresses are necessary—one 
for the internal memory and the other for the drum. Hence a complete 
word will be used for a drum order. A complete control cycle consists of 
3 parts: (1) an instruction cycle, where a (next) pair of instructions is brought 
from the memory to the instruction register, (2) an execution cycle A, where 
the left hand instruction is executed and (3) an execution cycle B, where the 
right hand instruction is executed. 


Since the number of orders is much less than 2°, the complications of 
complete decoding and recoding are avoided and static control is obtained 
by partial decoding of the order part. 


There are four sequencing units similar to the one used in ORDVAC 
type machines. In each sequencing unit, 2 flip-flops control the micro opera- 
tions. Four different combinations of these 2 flip-flops are possible, viz., 
11, 01, 00, 10 and the micro operations referring to these combinations for a 
particular sequencing unit, the Add-Shift unit (Fig. 5) will be described here 
in detail as a typical example of control logic design. Table I shows the 
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logic and polarity of the signals coming from the static control to this Add- 
Shift sequencing unit. 








AS COMP SIGNAL 















































































































age. es () 
ES(4) 
. 10 on CARRY RELEASE 
ain >] ; é 
: RO) 2 
tect p12} Yessas eats 
has " | 
; } | A ct MQ UR 
ine | / fe (2) 
of | J ogre se Met oe” E80 ve my 
nove tie Re). Oivide ) 
ow | | 7" ~ RANK ppg ps - 
| | u 
1e€r | TE MG ALR 
ns. } o ES@) : Toma ue 
| STOP as R15) RC) ES(!) R (5) 
CYCLING IF COMP. 
00 CL MQ Le 
, &s (2) R (7) 
he & R«s) 
Gr) ce ACC LR icc eniia ci 
Co omtis 1 care Teen ae — =—_ TOACC LR OT 
00 TO FALSE RANK R(8) 
Te ACC ur 
apt p — ty —~ ee, 
1 q 
8 OL—eq— GATE ACC uk 
J TO ACT LA at 
y GATa me uP Rim) 
» GaTe m@ UR Tome LRT 
/Y TO MQ.LR RT 
. Fic. 5. Add-Shift Sequence A-S. 
ie 
€ A-S sequence starts with 10 when it receives a signal from the previous 
f sequencing unit, (E.S. 4) it signals the Carry Release Driver. Here it may be 





mentioned that in order to have fast dying down of the carry level and to 
avoid persistence of the end-around-carry level in the case of division, carry 
to all the stages of the adder is pulled down by means of a driver and is released 
as A-S moves and is again pulled down when gating from Adder to accumu- 
lator Upper Register is complete. The same signal changes the state of the 
sequencing unit to 11 after a delay as shown in the ,diagram. Depending 
on R,, carry delay is initiated and 11 gives signal only when this delay is over. 
This signal goes to Clear Accumulator Upper Register Driver and clear M.Q. 
U.R. driver and after a delay changes the state of the sequencing unit to 01. 
This delay determines the pulse width of the clear pulse which in this computer 


has been made approximately 2p sec. 
















The gate pulse, depending on Rg, R; and R,, also starts at the same time 
as the clear pulse but is made to stay on for nearly ly sec. longer than the 
latter due to the delay shown after the 11 ‘ And” gate. 
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The 11 ‘ And’ gate feeds the pulse to the shift counter to keep track of the 
number of shifts. This shift counter* is a double rank counter similar to 
the one used in ORDVAC but incorporating an improved gating logic. A 
‘* recognition ” circuit associated with the shift counter, gives the “‘ recogni- 
tion”’ signal to signify the end of the Add-Shift operation. Thus in the 
case of shift orders, the recognition is obtained as soon as the number of 
shifts, as specified in the address part of the instruction, has been accom- 
plished and after 40 shifts in the case of multiplication and division orders. 
In the case of Add, Subtract and Collate orders the A-S sequence is allowed 
to go through its cycle only once. 


As in the case of 11, the clear Acc, L.R. and M.Q. L.R. and the corres- 
ponding gate pulses are obtained from the state 00 of this sequencing unit. 
The complete micro operations of the A-S sequencing unit can be listed as 
follows :— 








The state of A-S Unit Micro operations performed 
1 O1 .. (a) release the carry pulling. 
(5) delay the change of A-S to 11 if R, present (if 
R, —ve). 
1 | .. (a) clear Acc. U.R. 


(b) clear M.Q. U.R. 


(c) if Ry, gate adder output to Acc. U.R. In case 
of divide and no end around carry, write | 
in 2-4° (M.Q.) U.R. 


(d) if R;, gate Acc. L.R. to Acc. U.R. 

(e) if Rs, gate M.Q. L.R. to M.Q. U.R. 

(f) gate shift counter from False rank to True rank. 
0 1/0 .. (a) apply the carry pulling. 


(b) goes to the other sequencing unit (Execution 
Sequence) where it is used to stop A-S cycling 
depending on recognition. 


0 0 .. (a) Cl Acc. LR. 
(b) Cl. M.Q. L.R. if Rg. 
(c) if Rg, gate Acc. U.R. to Acc. L.R. (Down). 
(d) if Rg, gate Acc. U.R. to Acc. L.R. (Right). 
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he 
to § The state of A-S Unit Micro operations performed 
ee 
“A : (e) if Ry, gate Acc. U.R. to Acc. L.R. (Left). 
of &f (f) if Ry, gate M.Q. U.R. to M.Q. L.R. (Right). 
n- ff (g) if Ry, gate M.Q. U.R. to M.Q. LR. (Left). 
is. & . 
ad (h) gate the shift counter from true rank to false rank. 
Taking the case of the multiply order with 2-*° 
bit equal to 0 (i.e., the sense bit = 0). 
> A-S 1 0/1 .. releases the carry pulling. 
t. It does not wait for carry delay. 
. f A-S 1 1 .. Clears Acc. U.R., clears M.Q. U.R. 
. j As R; is present from the logic in Table I, Acc. L.R. 
q is gated to Acc. U.R. and as R, is present M.Q. 
L.R. is gated to M.Q. UR. 
, A-S 0 0 .. Clears Acc. L.R., clears M.Q. L.R. because R, is 
d present. As Ry, and R,, are present, Acc. U.R. 
f ‘ and M.Q. U.R. are gated “Right-Down” to Acc. 
: L.R. and M.Q. L.R. respectively. 
So the whole A-S cycle shifts the contents of Acc. 





and M.Q. Lower Register one place to the right. 





Taking the case of multiply with 2- = 1 (i.e., sence bit = 1), it can 
be seen that as R, is present, the carry delay will be initiated and as R, is 
present in place of Rg, after the carry delay, A-S 11 will gate the contents of 
the adder to Acc. U.R. Other micro operations will be same as in multiply 
with 2- — 0. 


So the whole A-S cycle adds the contents of Memory Register to Acc. 
and then shifts the contents of Acc, and M.Q. Lower Register one place to 
the right. 








The other sequencing units perform their micro operations in an exactly 
similar fashion. In addition to the Add-Shift unit described above, there 
are 4 other sequencing units, one for the instruction cycle, one for the execution 
cycle, one for the output control and the last one called the Main Control, 
to control all the sequencing units, 
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CLEAR AND GATE DRIVERS 


As mentioned earlier, clearing of the flip-flops(Fig. 2)is done by lowering 
the plate voltage of the left hand tube, thus lowering the grid of the right hand 
tube to cause it to cut off, and the left hand tube to become conducting. The 
gating is done by raising the common cathode of the flip-flop to cause the 
left tube to cut off. When the cathode is again lowered, the volage of the 
right grid is higher than that of the left grid. The right tube now is conducting 
and the left one is cut off. Clear and gate driver circuits‘ are similar to those 
described by J.C. Chu of Argonne National Laboratory for the Oak Ridge 
Automatic Computer. In order to minimise fluctuations of the driver output 
voltages, self-regulating drivers are used. 


For True/Complement gate and for carry pulling, drivers (+ 40 V)/ 
(—30 V) are used. They are also of the self-regulating type. Depending 
on the input, the output is constant either at + 40 V or — 30 V. Figure 6 
shows the circuit diagram of this driver. 
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Fic. 6. +40/—30 Driver. 
THE MEMORY 


The internal memory of the T.I.F.R. computer consists of a 3-dimensional 
ferrite core matrix of 1024 words each of 40 digits. Figure 7 shows the sche- 
matic diagram of the memory. 


As seen from the diagram, the memory consists of—apart from the store 
— X and Y diode decoding matrix for address selection, 40 output ampli- 
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fiers with amplitude and time discrimination for gating the output of memory 
cores into the Memory Register of the Arithmetic Unit, and read and write 
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Fic. 7. Schematic Diagram of Memory. 


circuits. The input to the memory can be either from the accumulator or 
the M.Q. or the Memory Register. In the case of the Acc, either the contents 
of the complete register or only the address part (either the left hand or the 
right hand address alone or both) can be written in. Since the memory read 
out is always destructive, in case it is required to rewrite the contents back, 
this is done from the Memory Register after the contents have been gated 
into it. 

Read and write logic is essentially the same as that described by Papian 
for the M.I.T. computer.’ Figure 8 shows the X-read, Y-read, Strobe, Write 
and Z (Inhibition) pulses. As can be seen the Y-read pulse starts after a 
slight delay from the X-read and the strobing is done after a further delay. 
The timings are shown in the figure. 


Figure 9 shows one of the four constant current read and write Drivers 
along with its pulse transformer. The diode matrix decoder supplies + 100 V 
at the selected terminal and + 60 V at the unselected terminals. These 
potentials are fed to the grids of the driver valves having pulse transformer 
in the plate circuit. All the driver cathodes through small potentiometers 
are joined together and kept at + 130 V ty a feed back stabilised system. 
All the driver valves are therefore normally cut off and no current flows 
through the pulse transformer. A current, slightly more than the driving 
current required, flows through the constant current tube. The read or write 
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pulse of about 30 V at the input, pulls down the grid of the constant current 
tube to 60 V. The common cathode is, however, arrested by the selected 
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Fic. 9. Memory X and Y Read Write Drivers. 


driver and the driving current flows through the selected pulse transformer. 
The magnitude of the driving current is almost independent of the variation 
in the driving tube characteristics. Minor variations in the currents through 
the individual transformers can be obtained by means of the individual cathode 


potentiometers of 100 ohms shown in the figure. 


The pulse transformers are wound on C cores obtained from the English 
Electric Company, London. The primary turns for read and write are 
125 each. They are wound on each leg of the C core and the secondary turns 


of 20 is wound 10 on each leg in series, 
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Figure 10 shows the circuit diagram of the pulse amplifier for the memory 
output. The rectifier diodes are biased slightly positive through the potential 
division at the cathode of the second stage of the amplifier. This helps in 
getting better signal to noise ratio. After the strobing at the 6AS6 ‘ And’ 
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Fic. 10. Amplifier for Memory Output. 


gate, the signal is fed to a univibrator to get a pulse of 4 sec. for reliable 
gating to the memory register. 


SUMMARY 


The T.I.F.R. computer is a parallel, binary, asynchronous single address 
machine. The arithmetic unit consists of a Memory Register, two (double 
rank) shift registers, the Accumulator and the M.Q. Register respectively, 
an Adder and a pair of True/Complement gates. A carry bypass logic is 
described which reduces the maximum carry delay time by a factor of 4. 
The maximum carry delay time for 40 bits is 3 sec. The clear and gate 
drivers for the Registers and the driver (+ 40 V/—30 V) for True/Complement 
gating are of self-regulating type and their design and operation are dealt with 
in detail. The Add-Shift sequencing unit is described as a typical example of 
the control logic design. Other sequencing units perform their micro opera- 
tions in exactly similar fashion. Schematic diagram of the fast access ferrite 
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core memory (access time 15 sec.) and some of the associated circuits are” 
described. 
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